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ABSTRACT 

We derive the physical properties of 580 molecular clouds based on their 
12 CO and 13 CO line emission detected in the University of Massachusetts-Stony 
Brook (UMSB) and Galactic Ring surveys. We provide a range of values of 
the physical properties of molecular clouds, and find a power-law correlation 
between their radii and masses, suggesting that the fractal dimension of the ISM 
is around 2.36. This relation, M = (228±18) F| 2 - 36 ± ao4 ; allows us to derive masses 
for an additional 170 GRS molecular clouds not covered by the UMSB survey. 
We derive the Galactic surface mass density of molecular gas and examine its 
spatial variations throughout the Galaxy. We find that the azimuthally averaged 
Galactic surface density of molecular gas peaks between Galactocentric radii of 
4 and 5 kpc. Although the Perseus arm is not detected in molecular gas, the 
Galactic surface density of molecular gas is enhanced along the positions of the 
Scutum-Crux and Sagittarius arms. This may indicate that molecular clouds 
form in spiral arms and are disrupted in the inter-arm space. Last, we find 
that the CO excitation temperature of molecular clouds decreases away from 
the Galactic center, suggesting a possible decline in the star formation rate with 
Galactocentric radius. There is a marginally significant enhancement in the CO 
excitation temperature of molecular clouds at a Galactocentric radius of about 
6 kpc, which in the longitude range of the GRS corresponds to the Sagittarius 
arm. This temperature increase could be associated with massive star formation 
in the Sagittarius spiral arm. 
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Introduction 



In the past forty years, there has been a considerable effort to establish the physical 
properties of molecular clouds, to understand how turbulence, gravity and magnetic fields 
shape their complex structure, how stars are born from their collapse, and if/how they relate 
to the spiral structure of the Milky Way. Several questions however remain unanswered. In 
particular, it is not clear how the spiral structure of our Galaxy relates to the formation 
and distribution of molecular clouds. This lack of evidence for a spiral structure traced by 
molecular clouds comes from a confused view of our own Galaxy, and can be accounted for 
by several difficulties in mapping spiral tracers in the Milky Way. In particular, distances to 
molecular clouds have proven very challenging to determine due to the kinematic distance 



ambiguity (ISanders et al.lll985t ISolomon et al .11 198 71 : IClemens et al. 



motions affecting the accuracy of kinematic distances ( iReid et al. 



1988), and to streaming 
20091 ). In addition, the 



lack of sampling, resolution, and sensitivity of CO surveys ha s made the deriva tion of the 
physical properties and structure of molecular clouds difficult ( IHeyer et al.l 120091 ). 

The lifetimes of molecular clouds also constitute a subject open to debate: are they 
long-lived (lifetime > 10 8 years) or short-lived (lifetime of 10 6 - 10 7 years) ? We still do 
not know whether m olecular clouds are transient features or gravitationally bound struc- 
tures. Simulations by iDobbs et al.l ( 120061 ) suggest that molecular clouds are unbound. Only 
clumps embedded in molecular clouds are gravi tationally bound and are subsequently able to 
form stars. On the other hand, simulations by lTasker fc Tan! (120091 ) suggest that molecular 
clouds are gravitationally bound. The gravitational state of a molecular cloud is described 
by its virial parameter. The virial parameter essentially represents the ratio of kinetic to 
gravitational energy in a molecular cloud. If it is greater than 1, the molecular cloud is not 
gravitationally bound. If it is smaller than 1, the molecular cloud is gravitationally bound. 

The physical properties of molecular clouds (mass, size, density, temperature, virial pa- 
rameter...) should not only reflect their formation and dynamical evolution, but also provide 
key information about the spiral structure of the Milky Way. Theoretical models predict that 
warm, diffuse gas is compressed in spiral ar ms, leading to the formation of atomic cloud s that 
give rise to molecular clouds within 10 My ( IDobbs et al.ll2006l ; iGlover fc Mac Lowll2010l ). The 
onset of star formation in high density regions cr eated by turbulent d ensity fluctuations and 
shocks occurs within a few global free fall times ( jKlessen et al.l 120001 ) . By the time a molec- 
ular cloud has formed and star formation been triggered, the molecular cloud has already 
reached the leading edge of the spiral arm where it formed, and starts to enter the inter- 
arm space. In the meantime, massive star formation photo-evaporates, photo-ionizes, and 
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dynamically disrupts the molecular cloud. As a result, molecular clouds should be rapidly 
(within a few million years) disrupted in the inter-arm space. If these theoretical predictions 
are correct, molecular clouds located inside spiral arms should be somewhat gravitationally 
bound while inter-arm molecular clouds should exhibit signs of disruption due to massive 
star-formation activity and dynamical disruption. 



Physical properties of molecular clouds have previously been investigated by lSolomon et al. 

f)l987h using the University of Massachusetts-Stony Brook ( UMSB) 12 CO J = 1 -> Galac- 
tic Plane Survey ( panders et al.l Il986t IClemens et al.lll986t ). Being undersampled in solid 
angle (with a resolution of 48" and a grid spacing of 3'), the UMSB survey did however 
not yield an accurate derivation of the physical properties of molecular clouds. In addition, 
molecular clouds detected in the UMSB survey suffered from blending, particularly near the 
tangent point, due to the use of an optically thick tracer ( 12 CO ). The Boston University 
- Five College Rad io Astronomy Observatory (BU-FCRAO) Galactic Ring Survey (GRS; 
Jackson et al.ll2006l ) is the first large scale 13 CO J = 1 — > survey that is both fully sampled 



and has an angular resolution sufficient to resolve the small scale structure of molecular 
clouds (resolution of 0.2 pc at a distance of 1 kpc). With a resolution of 46", a grid spacing 
of 22", a sensitivity of 0.13 K, and a coverage of 75 square degrees, it is a reliable data set 
to probe both the structure of molecular clouds and the Galactic distribution of molecular 
clouds. In addition, because the GRS uses an optically thin tracer (r ( 13 CO) ~ 1), and offers 
a spectral resolution superior to previous surveys, it does not suffer from blending as much 
as the UMSB, and can be used to reliably identify molecular clouds. 

In this paper, we derive the physical properties (mass, radius, surface density, ...) and 
the Galactic mass dist ribution of a sample of 750 molecular clouds id entified in the GRS 



( Rathborne et al.ll2009l ). for which kinematic distances are available from lRoman-Duval et al. 



(120091 ) . The paper is organized as follows. The data, methodology and conventions used to 
derive the physical properties of 580 molecular clouds detected in the GRS and covered by 
the UMSB survey are described in Sections [2] and [21 and the catalog is shown in Section 13.91 
In section HI we discuss the fractal structure the ISM, which can be derived from the relation 
between their sizes and masses. This relation allows us to derive masses for an additional 
170 GRS molecular clouds not covered by the UMSB survey in Section [SJ The histograms 
and range of values for the physical properties of molecular clouds are derived in Section 
El In Section [7J we derive the Galactic surface mass density of molecular gas and how it 
correlates with different models of the spiral structure of the Milky Way. Section [8] examines 
the variations of the physical properties of the GRS molecular clouds with Galactocentric 
radius. Finally, the limitations and caveats of this analysis are discussed in section [91 Section 
[TU] concludes this analysis. 
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Molecular line data 



2.1. 



13 



CO J = 1 -»■ GRS Data 



The sample of molecular clouds used i n this analysis was i dentified by their 13 CO J = 
1—7-0 emission in the BU-FCRAO GRS (I Jackson et al.ll2006l ). The GRS was conducted 
using the FCRAO 14 m telescope in New Salem, Massachussetts between 1998 and 2005. 
The survey, which used the SEQUOIA multipixel array, covers the range of Galactic longi- 
tude 18° < £ < 55.7° and Galactic latitude —1° < b < 1°. The survey achieved a spatial 
resolution of 46", sampled on a 22" grid, and a spectral resolution of 0.212 km s -1 for a 
noise variance a (T^) = 0.13 K [pT mb = 0.26 K accounting for the main beam efficiency of 
0.48). The survey covers the range of velocity —5 to 135 km s _1 for Galactic longitudes £ < 
40° and —5 to 85 km s" 1 for Galactic longitudes £ > 40°. 

Because the GRS uses 13 CO, which has an optical depth 50 times lower than that of 
12 CO , it allows for a better detection and separation of molecular clouds bot h spatially and 



spect rally than previous CO surveys. Using the algorithm CLUMPFIND ( 1 Williams et al. 



19941 ) applied to the GRS data smoot hed to 0.1° spatially a nd to 0.6 km s _1 spectrally, 
829 molecular clouds were identified by lRathborne et al.l (120091 ) . CLUMPFIND identifies as 
molecular clouds a set of contiguous voxels (i.e., (£, b, v) positions) with intensity values 
higher than a given threshold, which has to be determined empirica lly so as to best identify 
molecular clouds at all levels of emission. We refer the reader to iRathborne et al.l (120091 ) 
for the details of the identification procedure. Molecular cloud parameters such as Galactic 
longitude, latitude, and velocity of the 13 CO emission peak were estimated by CLUMPFIND. 
Individual 13 CO data cubes extracted from the GRS, and covering only the Galactic longi- 
tude, latitude, a n d velo city range of each individual molecular cloud were also created by 



Rathborne et al.l ( 20091). We use kinematic d istances t o 750 out of the 829 GRS molecular 
clouds derived in iRoman-Duval et al.l (120091 ) using the IClemensI ( 119851 ) rotation curve and 



H I self-absortion to resolve the kinematic distance ambiguity. 



2.2. 



12 



CO J = 1 UMSB data 



We make use of the 12 CO UMSB survey (jSanders et al.lll986l ; IClemens et al.lll9861 ). a 
joint program between FCRAO and the State University of New York at Stony Brook per- 
formed between November 1981 and March 1984. All of the observations were obtained using 
the FCRAO 14 m telescope. A grid sampled every 3' covering the range 18° < £ < 55° and 
— 1° < b < +1° was observed with a velocity resolution of 1 km s -1 and a spatial resolution 
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of 44". The sensitivity of the observations is 0.4 K per velocity channel. The UMSB survey 
only covers the velocity range -10 km s _1 < Vlsr < 90 km s^ 1 , which excludes the GRS 
molecular clouds located at low longitude and at the tangent point. 



3. Derivation of the physical properties of GRS molecular clouds 

3.1. Excitation temperatures 

The 12 CO J = 1 — > excitation temperature of our sample of 580 molecular clouds was 
derived from the 12 CO brightness temperature, with the assumption that the 12 CO line is 
optically thick (i.e r ( 12 CO) ^> 1). In the optically thick regime, the observed 12 CO brightnes s 



temperature T 12 and the excitation temperature T ex are related by iRohlfs &: Wilson! ( 12004 ): 



T K , = 5.53 ; r- (V, 



ln[l + 



5.53 



Ti2+0.837 



where both temperatures are in K. This equation includes a background subtraction ac- 
counting for the cosmic microwave background at T = 2.73 K. The excitation temperature 
was calculated at each voxel (i.e., (£, b, v) position) associated with each molecular cloud in 
the UMSB data resampled to the GRS grid. The 12 CO excitation temperature was used as 
a proxy for the 13 CO excitation temperature. The 13 CO and 12 CO excitation temperature 
should be identical in local thermal equilibrium since the energy levels of the two isotopomers 
are roughly the same. In reality, density gradients in molecular clouds and non-local ther- 
mal equilibrium situations can lead to differences in the excitation temperatures of 13 CO and 
12 CO. These effects are described in Section [9731 The mean excitation temperature of each 
molecular cloud was computed by averaging the excitation temperature over all the voxels 
that contain a 13 CO brightness temperature greater 4a or, accounting for the main beam 
efficiency, a minimum brightness temperature T13 > 1 K. 



3.2. 13 CO optical depths 



The 13 CO optical depth was derived using i) the 12 CO excitation temperature as a 
proxy for the 13 CO excitation temperature, and ii) the GRS 13 CO data. Since the UMSB 
and GRS data have nearly the same beam-width, corrections to account for different beam 
dilution between the tw o surveys can be neglected. The 13 CO optical depth is thus given by 



Rohlfs fc Wilson! (J200J): 
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Ti 3 (£, b, v) = -In 1 



0.189 T 13 (£,b, v) 



(2) 



0.16 



where T13 is the 13 CO optical depth and T13 is the background subtracted 13 CO brightness 
temperature in K. The 13 CO optical depth is evaluated at each GRS voxel (£, b, v) associated 
with each molecular cloud (i.e., in the individual data cubes associated with each molecular 
cloud). The line center optical depth was computed at each pixel (i.e., position (£, b) on the 
sky) where the 13 CO integrated intensity is greater than 4cr = §.23^/1% K km s , where N v 
is the number of channel in the cube. We also computed the average of the line center optical 
depth over all GRS pixels associated with each cloud where the 13 CO integrated intensity is 
greater than 4cr. 



3.3. 13 CO column densities 



From the optical depths and excitation temperatures, 13 CO co 



umn densities were de- 



rived at each GRS pixel (£, b) associated with each molecular cloud (IRohlfs fc Wilsonll2004j ): 

N( 13 CO)(e,b) ^ 26; , 1Q i4 f T ex {l,b,v)rM dv 
cm" 2 ' J 1 _ & T ^(d,v) kms^ 1 

13 CO column densities are computed by integrating Equation |3] only over GRS voxels as- 
sociated with each cloud where the 13 CO brightness temperature is greater than 4<r = 1 
K. 



3.4. Radii of molecular clouds 

The solid angle subtended by each GRS molecular cloud is computed by counting the 
number N pix of GRS pixels (i.e., (£, b) positions on the sky) associated with each molecular 
cloud that contain a 13 CO integrated intensity greater than 4cr = 0.23a//V^ K km s _1 : 

n = N pix AlAb (4) 

where A£ and Ab are the angular width of the pixels. Knowing the distance d to the molecular 
clouds, the area A and linear equivalent radius R of the molecular clouds are given by: 



A = ttd 2 



(5) 



3.5. Velocity dispersion 



The one-dimensional, intensity weighted, velocity dispersion of a molecular cloud is 
defined as: 

2 T,T 13 (eM>iK T M b ' v )( v - < v >) 2 m 

where T13 is the 13 CO brightness temperature in K, and only positions and velocity channels 
where the 13 CO main beam temperature is greater than 4<r = 1 K are taken into account. 
The average radial velocity < v > is given by: 

Er 13 (i,M>i/f T i3^M^ 
<v >= ^= rp To 7 — r (°J 



3.6. Masses 



Molecular hydrogen (H 2 ) and Helium (He) are the main constituents of molecular clouds. 
In order to derive molecular cloud masses from CO observations, one must therefore make 
assumptions about the abundance of CO relative to H2. First, we assume that the abun- 
dance of 13 C0 relative to H 2 and He is uniform. In reality, the CO abundance declines 
steep ly with decreasing Ay du e to photo-dissociation by the Galactic radiation field at 
< 3 (IGlover fc Mac Lowll2010l ). In constrast, self-shielded H 2 can exist at Ay as low as 
0.2 ( Wolfire et al.l 2010 ). As a result, the 13 CO /H 2 abundance is likely to be lower in the 
envelopes of molecular clouds than in their denser cores. The effects of abundance variations 
on the estimation of the molecular cloud masses are discussed in Section 19.21 Nevertheless, 
for simpli city, constant ratios n('^ 2 CO)/n(' 13 CO) = 45 and n( 12 CO)/n(H 2 ) = 8 x 1CT 5 are 



assumed (ILanger &: PenziasHl990l : iBlake et al.l 119871). A mean molecular we ight of 2.72 ac- 
counts for the presence of both H 2 and He (j Allen! 1 19731 : ISimon et al.ll200ll ). Under those 
assumptions, molecular cloud mas s es are derived from equation [9] below using kinematic 
distances from iRoman-Duval et al.l (120091 ) . 
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In this equation, the integrand is summed over all (£, b) positions where the 13 CO inte- 
grated intensity is greater than 4er. Along those sightlines, only velocity channels where the 
brightness temperature is greater than 4<r = 1 K are included in the integration. The second 
criterion is necessary due to numerical issues with Equation [TJ |2j and M when the brightness 
temperature value is negative. Hence, we do not include contributions from the noise in 
the integration described by Equation HI However, some noise peaks in voxels isolated from 
molecular cloud emission, with values greater than 4a, still remain even after filtering the 
integration with the first criterion. To remedy this problem, we filter the line-of-sights (£, b) 
along which the integration is perform by applying a threshold in 13 CO integrated intensity 
(first criterion). 



3.7. Number density and surface mass density 

The mean number density of particles (H2 and He) in the molecular clouds was estimated 
assuming spherical symmetry via: 

n(H 2 + He) M {4 R'y 1 

15.1 X — — X -7T — - (10) 



cm 3 M & \3 pc 3 

The surface mass density £ of the molecular clouds is defined as: 



M \ ( A x 1 



M pc~ 2 \M Q J Vpc 2 , 
where the area A of the molecular clouds is calculated as described in Section 13.41 



3.8. Virial mass and virial parameter 



Molecular clouds are supported against gravitational c ollapse by various mechanisms, 



such as turbulence, thermal gas pressure, and magnetic fields ( jHeitsch et al.ll200ll : lKlessen et al. 



20001 ) . Observations show that the linewidths of molecular clouds are much wider than their 



thermal line-widths. Transonic/supersonic turbulenc e must therefore be the main sourc e of 
kinetic energy and support in molecular clouds (e.g., lLarsonlll98ll : IWilliams et al.ll2000h . 

The virial parameter a of a molecular cloud is the ratio of its virial mass M„j r to its 
mass. It describes the ratio of internal, supporting energy to its gravitational energy. The 
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virial mass of a molecular cloud is defined as the mass for which a molecular cloud is in virial 
equilibrium, i.e., when the internal kinetic energy K equals half the gravitational energy U 
(2K + U = 0). It is given by: 

M mr = 1.3^ = 905- ^ (12) 
G pc (kms" 1 ) 2 

where R is the equivalent radius of a molecular cloud defined in section 13.41 and a VlD is the 
one-dimensional velocity dispersion defined in section [331 Note that the 3D isotropic velocity 
dispersion, a V) is v^3 times the ID velocity dispersion, which is measured in our spectroscopic 
data. For M > M vir (a < 1), 2K + U < and the molecular cloud is gravitationally bound. 
For M < M vir (a > 1), 2K + U > and the molecular cloud is not gravitationally bound. 

The virial mass is proportional to the linewidth and to the radius of a molecular cloud, 
with the proportionality constant depending on the number density profile. We have com- 
puted the CO number density profile by assuming spherical symmetry and deprojecting the 
column density with an Abel transform. By fitting a power-law to each den sity profile, we 



find that the average slope is —1.8. To be consistent with past literature ( ISolomon et al. 



19871 ). we use the proportionality constant derive from a density profile of slope —2 to com- 



pute the virial mass. 



3.9. Physical Properties of the GRS molecular clouds 

The first 25 entries of the derived physical properties of the 580 GRS molecular clouds 
covered by the UM SB are listed in Table [TJ The complete Table (including molecular clouds 
not covered by the UMSB survey, see Section EJ), can be found online. The first four columns 
indicate the molecular cloud name, Galactic longitude, latitude and LSR velocity from the 



Rathborne et al.l (120091 ) catalog. Columns 5 indicates the FWHM velocity dispersion of the 
molecular clouds (derived from section [3751 via a conversion factor of y/9>ln(2) between the 1- 
a and FWHM velocity dispersion). Column 6 gives the physical radii as defined in section 
[3j Columns 7 and 8 indicate the masses of the molecular clouds and their uncertainty. The 
derivation of error bars on the mass estimates are discussed in detail in Section 19.41 Column 
9 provides the mean number density of the molecular clouds. Columns 10 and 11 provide 
the mean excitation temperature and the mean 13 CO center-of-line optical depth. Column 
12 and 13 show the mean surface density of each cloud and the virial parameter. Column 
14 contains a flag "i" indicating that the cloud is covered by the UMSB survey. 
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On the fractal dimension of molecular clouds 



A dew decades ago , molecular clouds were thought t o be isolated, well-defined objects 
formed by coalescence (jOortl Il954j ; iField &: Saslawl Il965l ) and sustained in equilibrium by 
pressure from the hot inter-cloud medium. Since the 1980s, it has been known that molec- 
ular clouds are in fact dense sub-structures (n(Hg) > a few hundred cm" 3 ) in an underly - 



ing turbulent, fractal multi gas phase ISM (jScalolll985l . Il988c lFalgarondll989l ; IScalol Il990h . 



Molecular clouds form in gas over-densities resulting from supersonic turbulent flows. When 
these over-densities reach a visual extinction greater than a few tenths, molecules, such as 



li? and CO, become shielded from the p 


loto-dis 


sociating interstellar radiation field by dust 
[n this respect, molecular cloud boundaries 


(Glover & Mac Low 


2010; 


Wolfire et al. 


2010L 


are observed to be fractal ( 


Beech 


1987 


Scalo 


1990: 


Zimmermann & Stutzki 


1992). The 



interiors of molecular clouds ar e also fractal in n ature, as shown by their power- law den- 
sity spectra (IBrunt et all 12010), energy spectra (IRoman-Duval et al-llin preparation! ), and 
size-linewidth relation ( Larson! 1981 ; Heyer et al. 2009 ). On the smallest scales, very dense 
(n(H 2 ) > 10 5 cm -3 ) molecular cores give birth to star clusters. 

The fractal dimension of turbulent gas describes how completely it fills space as one 
zooms down to smaller and smaller scales. In other words, the fractal dimension corre- 
sponds to the degree of "sponginess" . The fract al dimension of molecular clouds has so far 



been i nvestigated via the perimeter-area relation (lBeechlll9871 ; IWakkerlll990l ; iBazell fc Desert 



19881 ; IScalol Il990l ; iFederrath et al.l 120091 ) . which relates the perime ter of a molecular cloud 
to its projected area on the sky. It has however been shown by iMandelbrot fc Whitrow 
(119831 ) that the radii R and masses M of sub-structures in a fractal are related via M oc 
R D . Because molecular clouds are sub-structures in an underlying fractal ISM, one should 
in principle be able to estimate the fractal dimension of the ISM within the range of spatial 
scales covered by molecular clouds from the correlation between the radii and masses of 
molecular clouds. Altough there are, in reality, other non-random mechanisms (e.g., spiral 
density waves) that modulate the distribution of cloud masses, treating molecular clouds as 
sub-structures in a fractal is a reasonable th order approximation. 

Figure [T] shows the correlation between molecular clouds' radii and masses. The radii 
and masses of GRS molecular clouds are related by a tight power-law correlation: M = 
228±18i? 2 36±0 04 , of exponent D = 2.36 ±0. 04. Within the error bars, this value is consistent 
with t he value of D derived i n the literature (IFalgarone fc Phillipslll991uElmegreen &; Falgarone 
19961 ; IFederrath et al.ll2009l ). The slope of the correlation was obtained by applying a chi- 
square minimization between a linear model and the observed relation between the logarithms 
of molecular clouds' radii and masses, weighted by the error on the mass. The error on the 
slope quoted here thus corresponds to the error on the linear fit between log(L) and log(M), 
including the error on the mass estimation for each molecular cloud. The error cited here is 
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quite small (4%) compared to the error cited by lElmegreen fc Falgarone (1996) for instance 
(30%) . The difference is due to the method used to derived D. lElmegreen &: Falgarone 
(119961 ) derived D from several surveys, some of which do not contain nearly as many molec- 
ular clouds as our sample. Hence, the errors quoted for the fractal dimension derived from 
individual surveys with sm all samples are larger th a n our error estimation. In addition, the 
final value of D quoted by lElmegreen fc Falgarone! (119961 ) corresponds to the average value 
yielded by all surveys, and the error on this value reflects the dispersi on between the different 



survey s, reduced with different methods and calibrat ions. However, lElmegreen fc Falgarone 
(1996) find D = 2.38±0.09 for galactic clouds from ISolomon et al.l ( 19871 ) and 



Dame et al. 



(Il986l ). which is consistent with both our value of D and our error calculation. 



The value of D between 2 and 3 corresponds to a spongy medium, which fills space 
more than simple sheets. This seemingly universal value of D, also observed in atmospheric 
clouds, could result from the intrinsic structure of supersonic, intermittent turbu lent flows 
(Sreenivasan fc Meneveau 1986: Meneveau fc Sreenivasan 1990: Sreenivasan 1991). 



5. Derivation of the masses of molecular clouds outside the UMSB coverage 

The tight correlation between the radii and masses of molecular clouds derived in Section 
IHcan be used to compute the masses of molecular clouds located outside the UMSB coverage, 
knowing their radii from Section 13.41 



M = (228 ± 18)i? 2 - 36±a04 (13) 

The density and surface density of those molecular clouds can then be derived using the 
method described in Section [3j We thus extended our catalog of molecular cloud masses, 
densi ties, and surface den s ities t o the 750 objects for which kinematic distances are available 



from iRoman-Duval et al.l (120091 ) . The online table includes GRS molecular clouds located 
outside the UMSB survey coverage, which are flagged by an "o" in the last column of the 
online table and have their temperature and optical depth set to zero. Molecular clouds 
covered by the UMSB are flagged by an "i" . Figure [2] shows the spatial distribution of GRS 
molecular clouds covered by the UMSB survey, and outside its coverage. 



6. Histograms of molecular clouds' physical properties 



This section is intended to provide a range of values for the physical properties of molec- 
ular clouds, that may later be used in other studies. For instance, they may be important 
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to constrain models of molecular cloud formation and evolution. The masses and radii of 
molecular clouds might also be used to predict g amma ray fluxes emanating from the inter - 
action between molecular clouds and cosmic rays (lAharonian et al.ll2008t iGabici et al.ll2009l ). 
Furthermore, the comparison of the mass spectrum of molecular clouds to that of clumps and 
cores is essential to understand the fragmentation process that leads from molecular clouds 
to stars. Last, the radius and mass s pectra of molecular clouds hav e been shown to result 
from the fractal structure of the ISM (jElmegreen fc Falgarondll996l ). They are therefore of 
great interest if one wants to constrain the fractal dimension of the ISM. The histograms of 
the physical properties (mass, radius, density, velocity dispersion, surface mass density, and 
virial parameter) of the GRS molecular clouds are shown in Figures [3] and HJ 



6.1. Radius and mass distributions 



The top panels of Figure |3] show the radius and mass spectra of the sample of 750 GRS 
molecular clouds, ^(R) = dR/dln(R) and $(M) = dN/dln{M). For M > 10 5 M , the mass 
spectrum follows a power law: $(M) oc jVf- 1 - 6 ^ - 25 . The radius spectrum also follows a 
power-law for R > 10 pc: ^(R) oc fj; -3 ' 90±0 ' 65 . Within the error bars, the slope of the mass 
spect rum derived in this paper is consistent with the value of —1.5 obtain ed in previous work 



e.g., 



Sanders et al. 1985 



Solomon et all 119871 : IWilliams fc McKedll997h . The slope of the 
radius spectrum is higher than the value obtained bv lSanders et al.l (1985T) . but is consistent, 



within the error bars, with the radius distribution obtained by lHeyer et al.l (120011 ). 

Our sample of molecular clouds is complete above the turn-over mass of M to = 4xl0 4 
M Q , such that the slope of the mass spectrum of molecular clouds should not be affected 
by a lack of completenesss. Our sample of molecular clouds was identified in a version of 
the GRS smoothed to 6' spatially and to 0.6 km s -1 spectrally. GRS molecular clouds were 
det ected by CLUMPFIND as contiguous voxels of brightness temperature greater than 0.2 
K (IRathborne et al.l 120091 ) . iRathborne et al.l (120091 ) furthermore applied the condition that 
molecular cloud candidates detected by CLUMPFIND must contain at least 16 smoothed 
voxels in order to be identified as a molecular cloud. Assuming that the 13 CO line is optically 
thin, the minimum mass of a molecular cloud is given by: 



M d 2 
> 0.05 5- T e , 

Mq ~ kpC 2 



x 16T„ 



A£Ab AV 
~> r kms- ] 



(14) 



where d is the distance in kpc, T ex is the excitation temperature in K, and T m j„ = 0.2 K 
is the threshold brightness temperature (corrected for beam efficiency). Hence, with T ex = 
6.32 K (the average value observed in the GRS), A£ = Ab = 0.1°, and AV = 0.6 km s -1 , 
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the completeness limit is M min = 50 rf| , where dk pc is the distance in kpc. Thus, M min = 
200 M Q at 2 kpc, 1250 M Q at 5 kpc, 5000 M at 10 kpc, and 11250 M & at 15 kpc. Since the 
maximum distance probed by the GRS is 15 kpc (based on the GRS Galactic longitude and 
velocity coverage), the turn-over mass is greater than the completeness limit of the GRS, 
and the slope of the mass spectrum of molecular clouds above the turn-over mass should not 
be affected by completeness effects. 

This calculation of the completeness limit does not take into account confusion, which 
in reality has an effect on our data. Confusion is more pronounced near the tangent point, 
where large physical separations correspond to small radial velocity differences. There is also 
the problem of molecular clouds that have similar radial velocities, but are located on either 
side of the tangent point, at the near and far kinematic distances. In these cases, if the 
velocity difference between the two clouds is small enough (typically less then the line-width 
of the clouds), the clump-finding algorithm used to identify the GRS clouds, CLUMPFIND, 
will blend these two clouds into a single object. The extent of the calculation and modeling 
required to estimate how confusion affects the completeness limit by far exceeds the scope of 
this paper. The blending of two molecular clouds depends on many parameters such as the 
geometry (their location in the Galaxy), the structure of large scale galactic features, the 
distance and radial velocity of the clouds (related by the rotation curve), their line- widths, 
radii, angular separation, and the parameters used in CLUMPFIND. Nonetheless, confusion 
in the GRS, which uses 13 CO as a tracer, is not nearly as severe as in surveys using the 
optically thick 12 CO as a tracer for molecular gas. 

It is worth mentioning that the GRS has a limited field-of-view, only covering the 
Galactic latitude range —1° < b < 1°. This in principle imposes a limit on a maximum 
cloud's radius, R ma x, and mass, M^^. y , detectable by the GRS. We can use the relation 
between a cloud's mass and radius derived in Section H] to estimate this maximum mass, 
with R max = vrrf/180. Thus, we find M™/ = 1.96 x 10 5 d 2 k ™. We will however show in 
Section I9TT1 that this limit is never reached, and thus does not affect our sample of molecular 
clouds and their physical properties. 

6.2. Virial parameter 

The fourth panel of Figure [3] shows the histogram of the 750 molecular clouds' virial 
parameter, the median value of which is 0.46±0.07. 70% of our molecular cloud sample (both 
in mass and number) have virial parameters < 1. This analysis thus suggests that most of 
the molecular mass contained in identifiable molecular clouds is located in gravitationally 
bound structures. 
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6.3. Number density and surface mass density 



The bottom left panel shows the mean density of H2 in our sample of 750 molecular 
clouds, the median of which is 231 cm -3 . This value is well below the critical density of the 
CO J = 1 — > transition, n cr = 2.7xl0 3 cm -3 , suggesting that the gas with density n > 
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n,. 



is not resolved by a 48" beam (0.25 pc at d — 1 kpc), and that its filling factor is low. 
The bottom right panel shows the surface mass density of the molecular clouds, with 
a median of 144 M c .) pc" 2 . Using the Galactic gas-to-dust ratio < Nu/Ay > = 1.9 xlO 21 



cm 2 mag 1 (IWhittetl 120031 ) . this corresponds to a median visual extinction of 7 mag. This 



valu e is consisten t with the prediction from photo-ionization dominated star formation the- 
ory (lMcKedll989l ). A median su rface mass dens i ty of 140 M pc~ 2 is lower than the median 
value of 206 M pc -2 derived by lSolomon et all ( 119871 ) based on the vi rial masses of a s ample 
of molecular clouds identified in the 12 CO UMSB survey. Note that ISolomon et all (119871 ) 
originally found a median surface mass density of 170 M & pc -2 , assuming that the distance 
from the sun to the Galactic center is 10 kp c. Assuming a Ga lactocentric radius of 8.5 kpc 
for the sun, this value becomes 206 M pc -2 (IHeyer et al.ll2009f ). T he median surface d ensity 
derived here is also higher than the value of 42 M & pc -2 derived by IHeyer et all (120091 ). who 



re-examined the masses and surface mass densities of the ISolomon et all ( 19871) sampl e using 



the GRS and a method similar method to ours. Similar to our analysis, IHeyer et all (120091 ) 
estimated the excitation temperature from the 12 CO line emission and derived the mass and 
surface dens ity from 13 CO GRS m easurements and the exci tation temperature . 

For the ISolomon et all ( 119871 ) molecular cloud sample, IHeyer et al.l (120091 ) found a me- 
dian surfac e density of 42 M^pc" 2 using the area Al (the 1 K isophote of the 12 CO line) 
defined by ISolomon et al.l (119871 ) to compute masses and surface mass densities. However, 
computing surface mass densities within the half power 12 CO isophote (A2) yie lds a median 



surface mass density close to 200 M Q pc 2 (see Fig. 4 of IHeyer et a 



2009). It is thus 



likely t hat the discrepancy b etween the surface densities derived here, in IHeyer et al.l ( 120091 ). 
and in ISolomon et al.l ( 119871 ) be explained by the different methods and thresholds used to 
compute t he molecula r clouds' properties. The 12 CO line is about 4-5 times as bright as the 
13 CO line (Lisztl 2006h . A s a result, it is likely that a significant fraction of the area Al used 
by lSolomon et al.l ( 119871 ) to compute the molecular clouds' masses using the 12 CO line only 
by assuming that molecular clouds are virialize d does not exhibit 13 CO line emission above 
the detection threshold of the GRS. Since the IHeyer et al.l ( 120091 ) derivation of molecular 
clouds' masses is based on 13 CO as a tracer of molecular gas column density, this would re- 
sult in a dilution of the surface mass density derived from 13 CO , which woul d appear lower 
compa red to the surface density derived from the brighter 12 CO emissio n by ISolomon et al. 
( 1987 ). This effects likely contribu t es to the median surface density in Heyer et al. ( 2009h 



being lower than in ISolomon et al.l ( 119871 ) and in our analysis. Since our derivation of the 
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molecular clouds' mass, radius, and surface mass density is only performed within the 4cr 
contour of the 13 CO integrated intensity, it is likely biases towards the most opaque regions 
of molecular clouds, where the surface mass density is higher th an in the envelope s, thus 
raising the median value of the surface mass density compared to iHeyer et al.l ( 120091 ) . 



6.4. Excitation temperature and optical depth 



Figure H] shows the histograms of the excitation temperature and 13 CO optical depth of 
the 580 GRS molecular clouds covered by the UMSB survey. The mean excitation temper- 
ature is 6.32±0.04 K, and the mean optical depth of the 13 CO line is 1.46±0.02. Although 
the 13 CO line is less oqaque than its 12 CO counterpart, optical depth effects should be taken 
into account in the derivatio n of 13 CO column densities from 13 CO spectral line mapping. 

Rathborne et al.l (120091 ) previously derived the excitation temperature and 13 CO opti- 



cal depth of the GRS catalog of molecular clouds, and found a lower mean optical depth 
(0.13), and a higher mean excitation tempe r ature (8.8 K). The difference between the re- 
sults pres ented here and in IRathborne et al.l (120091 ) can be explained by the use of different 
methods. IRathborne et al.l ( 120091 ) first computed the mean 13 CO and 12 CO brightness tem- 
peratures each GRS molecular cloud, < T 13 > and < T 12 > before applying Equations [1] 
and [2] to obtain the "mean" excitation temperature and optical depth, T ex (< T 12 >) and 
t 13 (< T 13 >, < T ex >). In contrast, we first computed the excitation temperature T ex and 
the 13 CO optical depth ri 3 from the 13 CO and 12 CO brightness temperatures using Equa- 
tions [1] and [2] at each voxel associated with a particular molecular cloud. Only then did 
we average the excitation temperature over voxels with 13 CO brightness temperature above 
the noise threshold (4<r), and the line center optical depth over pixels with 13 CO integrated 
intensity above the 4cr noise level to obtain < T ex (T 12 ) > and < r 13 (T 13 , T ex ) >. Because 
Equations [Hand [2] are non linear with T12, Ti 3 and T ex , the difference between performing the 
average before or after applying the equation can be quite large. For instance, for GRSMC 



G053.59+00.04, we find < r 13 (T 13 , T„) > = 1.78, wh i 



e ri 3 (< T 13 >, < T ex >) = 0.24, close to 



the typical optical depth fo und by IRathborne et al.l (120091 ). In this respect, molecular cloud 
masses recently derived by Urquhart et al. ( 201oh likely under-estimate the masses of GRS 



molecular c louds, since they use th e (uniform) excitation temperatures an d optical depths 



derive d by IRathborne et al.l (120091 ). We find that the masses derived by lUrquhart et al. 
(120 101 ) are lower than masses derived here by a factor 2 to 3. 
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Galactic mass distribution of molecular clouds 



Roman-Duval et al.l (120091 ) derived the Galactic CO surface brightness, which roughly 



represents the Galactic surface density of molecular gas traced by CO. This analysis, however, 
did not account for excitation temperature variations and 13 CO optical depth effects. Due to 
variations in the heating rate with local environment, the excitation temperature is however 
likely to vary significantly with Galactocentric radius and local environment. The densest 
regions of molecular clouds exhibiting 13 CO optical depths > 1, non-linearities between 13 CO 
surface brightness and molecular gas masses should be taken into account. Molecular cloud 
masses derived in this paper therefore allow for a more accurate and rigorous derivation of 
the Galactic mass distribution of molecular clouds. We derived the surface mass density of 
molecular gas in the Milky Way by summing the masses of GRS molecular clouds over circular 
bins of radius 0.5 kpc, with a sampling of 0.12 kpc. The resulting galactic surface mass 
density of molecular gas is shown in Figure [5) The green, yellow, red, and blue lines represent 
the 3 kpc arm, the Scutum-Crux arm, the Sagittarius arm, and the Perseus arm from the 
four-arm model by IValled (119951 ) respectively, which is based on a compilation of different 
tracers from the literature (e.g., CO, Hn regions, magne tic fields, ele ctron density etc.). 
The locations of the Scutum-Crux and Perseus arms from I Vailed (1995) are coincident with 



a two-arm model from Benjamin et al.l (120051 ); iBenjaminl (120091 ); IChurchwell et al.l (120091 ) 
based on K-giants and M-Giants star counts from the Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE). The Sagittarius arm is not detected as an overdensity 
i n the old stellar popula tion, perhaps because the Sagittarius arm is only a gas compression 
dChurchwell et al.ll2009l ). 

The surface mass density of molecular gas deduced from the GRS appears to be enhanced 
along the Scutum-Crux and Sagittarius arms. This suggests that, despite the uncertainty 
in kinematic distances due to non-circular motions near spiral arms and the uncertainty 
in the rotation curve, molecular clouds are good tracers of the Scutum and Sagittarius 
arms. Although the Perseus arm has previously been det ected in several tracers (e.g., water 



mase r s, molecular gas, FIR o bservations, star counts, see iHeyer fc Terebeylll998l ; iReid et al. 



20091 ; IChurchwell et al.l 120091 ). the Perseus arm is not detected as a strong enhancement in 
the Galactic surface density of molecular gas. This could be due to several effects. First, 
molecular clouds located at 10-15 kpc, on the far side of the Galaxy, tend to have smaller 
angular sizes than molecular clouds located at closer distances. Because the velocities of near 
and far molecular clouds located at the same galactocentric radius are the same, it is possible 
that far molecular clouds be blended together with near molecular clouds and assigned to the 
near kinematic distance. Second, non-circular motions near the Perseus arm may cause the 
distance estimate to molecular clouds located in the Perseus arm to be inaccurate. Indeed, 
there are two large molecular complexes located between the Sagittarius and Perseus arms, 
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which could be associated with the Perseus arm if their distance is systematically under- 
estimated. Third, the completeness limit of the GRS at 10-14 kpc is 5xl0 3 -10 4 M & , such 
that only a handful of very massive molecular clouds can be detected at the distance of the 
Perseus arm. Last, the molecular content of the Milky Way in the Galactocentric radius range 
of the Perseus arm (R ga i = 7-8 kpc) is of order T, ga i = 1-2 M e kpc -2 , much lower (by at least 
an order of magnitude, see Section |8T2|) . than in the 5 kpc molecular ring. The combination 
of these last two effects — completeness effects and low molecular surface density at R ga i 
= 7-8 kpc — can be further quantified. Assuming a power-law mass spectrum of exponent 
— 1.6 (see Section [6] and also 



Sanders et al. 1985 



Solomon et al. 1987 



Williams fc McKee 



19971 ). a Galactic surface density of molecular gas T, ga i = 1-2 M kpc 2 , a completeness 



limit of 10 4 Mq at a distance 10-15 kpc, and a maximum mass of 10 6 M Q (the maximum 
mass detected in the GRS), we predict that there should be 4 molecular clouds per kpc 2 
above the completeness limit in the region of the Perseus arm. Indeed, this number would 
not be sufficient to resolve the Perseus arm. This prediction is also supported by the actual 
number surface density of molecular clouds observed in the GRS in the Perseus arm region. 
The fact that the measured surface number density of molecular clouds in the vicinity of 
the Perseus arm agrees with our simple prediction based on the completeness limit and the 
azimuthally averaged Galactic surface density of molecular gas suggests that completeness 
effects and the low Galactic molecular content at Galactocentric radii 7-8 kpc likely account 
for the non-detection of the Perseus arm in the GRS. However, we do find localized egions 
with 0-1 cloud kpc -2 . Combined with localized enhancements of the number and surface 
mass density in the inter-arm space between the Sagittarius and Perseus arms, this could 
suggest that distance uncertainties due to non-circular motion near spiral arms play a role 
in the non-detection of the Perseus arm. 

The theoretical implications of the enhancement of t he Galactic surface mass d ensity 
of molecular gas along spiral arms have been discussed in iRoman-Duval et al.l ( 120091 ): the 
confinement of molecular clouds to spiral arms suggests that molecular clouds must form 
in spiral arms via a combination of hydrodynamical processes due to the compression in 
the spiral arm, self-gravity, Parker instability, and orbit crowding, and be disrupted in the 
inter-arm space. 
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8. Variations of the physical properties of molecular clouds with 

Galactocentric radius 

8.1. CO excitation temperature versus Galactocentric radius: effects of 

star-formation ? 



8.1.1. Decline in CO excitation temperature with Galactocentric radius 

The equilibrium temperature of molecular clouds results from a balance between heat- 
ing and cooling. The heating of gas in molecular clouds is due to: i) electrons ejected from 
dust grains illuminated by t he Interstellar Rad iation Field (ISRF) due to the photoeletric ef- 
fect 



Bakes fc Tielendll994 ; IWolfire et al.1120031 ). ii) H 2 photo-dissociation 



1977 ). iii) collisions between Ga lactic Cosmic Rays ( GCRs) and molecules ([Goldsmith fc Langer 



Black fc Dalgarriol 



978 ), iv) UV pumping of H 2 (IBurton et al.lll990l ). and v) formation of H 2 on dust grains 



( iHollenbach fc McKedll989l ). In the densest, most shielded regions of molecular clouds that 



are not penetrated by the ISRF, heating by GCRs dominates. In the regions closer to the 
CO boundary, photoelectric effect likely plays an important role in the heating of the gas. 

Fine structure lines of metals, most importantly [C n] at 158 /nn, and radiative CO ro- 
tational transitions, are responsible for the cooling of gas in molecular clouds. The decrease 
in the abundance of metals away from the Galactic center effect ively decreases the cooling 
rate in molecular clouds as their galactocentric radius increases ( jOuireza et al.ll2006l ). 

Because the gas temperature is the result of the equilibrium between heating and cool- 
ing, the effects of the variations in the heating and cooling rates with galactocentric radius 
on the gas temperature are reflected in the variations of the gas temperature throughout the 
Galactic plane. The CO excitation temperature however depends on both the gas tempera- 
ture and the local number density. T ex can be significanlty lower than the gas temperature 
if the density is lower than the CO critical density, in which case there are not enough 
particules to collisionally excite the CO molecule to the J=l level. This effect is known as 
sub-thermal excitation. Nonetheless, the CO excitation temperature is likely close to the gas 
temperature in the dense regions of molecular clouds, which dominate the emission. Hence, 
the variations of the CO exciation temperature with Galactocentric radius should constrain 
the balance between gas heating and cooling. 

Figure [6] shows the CO excitation temperature versus R ga i- Only molecular clouds cov- 
ered by the UMSB survey, for which we could derive a CO excitation temperature and optical 
depth, are taken into account. The top panel shows the maximum and mean excitation tem- 
perature in each molecular cloud. Since the mean excitation temperature is likely lower 
than the gas temperature due to sub-thermal excitation effects, the maximum excitation 
temperature in each cloud, corresponding to the densest regions, probably reflects the gas 
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temperature more accurately. The bottom panel shows these temperatures averaged over 
0.3 kpc galactocentric radius bins. In both cases, the CO excitation temperature decreases 
smoothly with galactocentric radius, from 14 K at R ga i ~ 4 kpc down to 8 K at R ga i ~ 8 
kpc (for the maximum excitation temperature). 

This large-scale decline in the gas temperature away from the Galactic center indicates 
that the slight decline in metallicity (and hence cooling rate) away from the Galactic center 
is not large enough to overcome the decrease in the heating rate (due to a decrease in the 
strength of the ISRF and/or the GCR flux) with Galactocentric radius. We have modeled 
and quantified the contribution of variations in the GCR flux and the ISRF to the varia- 
tions of the gas temperat ure with Galact o centri c rad ius. The variations of the GCR flux 
have been investigated by iBloemen et al.l (119861 ) and IStrong et al.l (119881 ) . While the flux 
of cosmic ray nuclei n appears constant with Galactocentric radius, the flux of cosmic ray 
electrons, </> R , s h ows a gradient best described by 4> e (R ga i) = ^(-Ro^" ' 19 ^"^-*. Note that 
Bloemen et al.l ( 119861 ) originally cite a gradie nt of —0.16 kpc" 1 , assum ing R Q = 10 kpc. 
With the more recent value of R Q = 8.5 kpc (IKerr fc Lynden-Belllll986l ). this corresponds 
to a gradient of —0.19 kpc" 1 . The heating rate by GCR is then Tqcr = (fie(Rgai) A cr n(_£f 2 ) 
( jGoldsmith fc Langerlll978l ). where A cr = 17 - 26 eV is the energy deposited as heat as a 



resu lt of the ionization by GCR s. The GCR flux at R ga i = R& is in the range 1.5 - 3 xlO 17 
s" 1 (jGoldsmith &: Langerlll978l ). Hence, the heating rate by GCRs is given by: 



GCIi(Rgal 



5 x 10- 28 n{H 2 )e-°- 19{R - R ^ [ergs cm" 3 s" 1 ] 



(15) 



where we chose a value of e A cr = 5x10 28 ergs s 1 , whi ch best describes the excita tion 



temperature trend in our data, and is in the range cite d by iGoldsmith fc Langerl (119781 ). 

The heating rate by photo-electric effect is given in Bakes fc Tielensl ~( 1994 ). For neutral 
grains, the heating rate is given by T pe = 4.86xlO~ 26 G (R ga i)n H [ergs cm" 3 s" 1 ], in an 
unattenuated medium, where G is the strength of the ISRF in Habing units (Go varies 
with Galactocentric radius), and rin is the hydrogen density. In a molecular cloud of visual 
extinction Ay, the heating rate due to photo-eletric effect is therefore: 



pe 



4.86 x 10- 26 G (R ga i)e- 5AAv n H [ergs cm" V 



(16) 



In this equation, we assumed that CO was present for Ay > 1 ( jWolfire et al.l 120101 ) . and 
used the standard Milky Way ext inction curve to rel ate A 10 oo, the FUV extinction at 1000 

he st rength of the ISRF, G , can 



A, and Ay (A 1000 = 4-5 Ay, see iGordon et~aD l2003h. 
be observed via the dust temperature dBernard fc et al. 



emissivity index of the dust ((3 = 1.5-2, see iBoulanger et al. 



2008]): Gn oc Tj +/3 , where is the 



1996; 



Gordon et al.ll2010h . and 
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Td is the equilibrium temperature of large dust grains. In addition, ISodroski et al.l (119971 ) 
derived the dust equilibrium temperature versus Galactocentric radius. They found T d (R ga i) 
= 28 K - 6.8 K x (R gal /R e ). Hence, with the local IRSF having strength G (R Q ) = 1.7, 
the variations of Gq throughout the Galactic plane are given by: 



G (R ga i) = G (R Q ) (l.32 - 0.32^ 



(4+/3) 

(17) 



This expression for Go(R ga i) can then be used in Equation [16] to calculate the heating due 
to photo-electric effect function of R ga i- 

The cooling rate by a variety of fine structure atomic and molecular lines is given in 



Goldsmith fc Langerl (I1978I ). For n(H 2 ) = 300 cm , which is close to the median density 



(^rned _ 23Q cm -3 ) in our sample of molecular clouds, the cooling rate is A = 4.7xlO _27 T 1,6 
[ergs cm -3 s -1 ]. 

In thermal equilbrium, the gas temperature can then be determined by equating the 
cooling and heating rates: A = Y. The gas temperature resulting from the balance between 
heating by GCR and line emission cooling is shown as a blue line in Figure [6j For the H2 
number density, n(H 2 ), we use the median value of our sample (n^ d = 230 cm -3 ) in equation 
[TBI Our simple calculation reproduces the variations of the CO excitation temperature with 
Galactocentric radius within the error bars. In Figure EJ we also show the gas temperature 
obtained from the equilibrium between cooling and heating by photo-electric effect (dashed 
green line). We used Equations [TBI a nd [T71 with the med ian number density found in our 



molecular cloud sample, and /3 = 1.5 ([Gordon et al.ll2010l ). In this case, our prediction also 



reproduces the temperature trend within the error bars. 

Thus, both a decline in the GCR flux and in the strength of the ISRF can explain the 
decrease in CO excitation temperature derived in this analysis based on the GRS. It is likely 
that the GCR flux and the strength of the ISRF are correlated, because they both depend on 
the star formation rate. In this case, the decline in gas temperature with R ga i would reflect a 
decrease in the star formation rate with Galactocentric radius. Indeed, young massive stars 
produce most of the ionizing and dissociating radiation responsible for gas heating in the 
envelopes of molecular clouds. GCRs are thought to be produced in supernova remnants via 
Fermi acceleration. Since supernovae remnants are the products of massive star-formation, 
the GCR flux is therefore tied to massive star-formation. Because of the large mean free 
path of GCRs, however, it is not clear how the heating rate due to GCRs in molecular 
clouds relates to the local star formation rate. Nonetheless, a gradient of — 30%/i? Q in the 
star formation rate (see Equation ITT)) would explain the decline of the gas temperature with 
Galactocentric radius. 
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8.1.2. Enhancement at R ga i = 6.4 kpc 

The maximum CO excitation temperature in each cloud is enhanced by about 2 K at 
a galactocentric radius of 6.4 kpc. This enhancement is significant at the 3a level with 
respect to the uncertainty plotted in Figure [6j Over the longitude coverage of the GRS, 
a galactocentric radius of 6.4 kpc corresponds to the Sagittarius arm. An elevated gas 
temperature at 6.4 kpc may be related to massive star-formation occuring in the Sagittarius 
arm (both through the strength of the radiation field and the GCR flux). In Figure El an 
increase of 50% in the GCR flux (blue line) and in the strength of the ISRF (green line) at 
Rgai = 6.4 kpc reproduces the trend observed in the GRS. 

The idea that star formation locally increases the gas temperature inside mo lecular 



clouds can be further tested by comparing different data sets. I Anderson et al.l ( 120091 ) found 
an enhancement of H II regions at a galactocentric radius of ~ 6 kpc. In addition, the 15 
hottest molecular clouds located at galactocentric radii 6-7 kpc (which stand out in Figure 
If)]) all have confirmed star formation activity - PAH emission in the 8 fim band, and 24 fim 
emission tracing warm dust- in the GLIMPSE and MIPSGAL infrared surveys. Eight of 
these fifteen molecular clouds are associated with the star formation region W51 located 
at a galactic longitude of 49° - 50° and at a distance of 5.5 kpc, which places it near the 
Sagittarius arm. The seven others belong to different star formation regions and contain 
known masers and H II regions. To test the hypothesis that star formation activity increases 
the gas temperature of molecular clouds, we compare the excitation temperatures of "active" 
m olecular clouds, conta i ning H II regions, to the temperature of "quiescent" molecular clouds. 



In iRoman-Duval et al.l (120091 ) . a catalog of such "active" and "quiescent" molecular clouds 
was established based on the coincidence between the morphology of the 2 1 cm continuum 
and 13 CO emission, and the velocity of 13 CO and recombination lines from I Anderson et al. 



( 120091 ). This catalog can readily be used to compare the excitation temperatures of active 
and quiescent molecular clouds. The histograms of the CO excitation temperature of both 
samples are shown in Figure [71 On average, "active" molecular clouds have a slightly higher 
temperature, < T ex > = 6.96±0.22, compared to "quiescent" molecular clouds (< T ex > 
= 6.28±0.04). A Kolmogorov Smirnov test shows that the active and quiescent molecular 
cloud populations have a less than 1% chance of being from the same distribution. This 
difference of temperature between active and quiescent molecular clouds, significant at the 
3a level, supports the idea that star formation increases the gas temperature (and hence CO 
excitation temperature) of molecular clouds. 
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8.2. Mass and surface density 

The variations of the masses and surface mass densities of the GRS molecular clouds as 
a function of galactocentric radius should help understand and predict the variations in the 
star formation activity resulting from this reservoir of molecular gas throughout the Galactic 
disk. The azimuthally averaged Galactic surface mass density of molecular gas, S gal , was 
obtained by summing the molecular clouds' masses over 0.5 kpc galactocentric radius bins. 
Since the range of longitudes covered by the GRS data is limited, each radial bin covers only 
a limited range in azimuth, which varies with Galacotcentric radius. To convert the mass 
measurements to a more uniform measure of gas content, the total mass contained in a radial 
bin was therefore divided by the surface area covered by the GRS data within the radial 
bin. As a result, we obtained an azimuthally averaged surface mass density of molecular gas 
encompassed by GRS molecular clouds at each galactocentric radius. The left panel of Figure 
|8] shows E ga i as a function of galactocentric radius R ga i. S ga i peaks at R ga i = 4.5 kpc (X ga i = 
2.5 xlO 6 M kpc~ 2 ), and decreases steeply with galactocentric radius, down to 10 5 M kpc~ 2 
at R ga i = 7.5 kpc. 

The middle panel of Figure [8] shows the molecular clouds' masses averaged over 0.5 kpc 
galactocentric radius bins. Since each Galactocentric radius bin covers a range of distances, 
and since the completeness limit varies with distance, only molecular clouds with masses 
greater than the completeness limit at 15 kpc, or 1.1x10 M pc~ 2 , were taken into ac- 
count. The mean molecular cloud mass in each radial bin appears rather constant with R ga \ 
out to 5.5 kpc, at which point the mean molecular cloud mass starts to drop from 10 5 M Q 
at i2g a i = 5.5 kpc down to 5xl0 4 M at _R ga i = 7.5 kpc. 

The right panel of Figure [8] shows the molecular clouds' surface mass density, S c , av- 
eraged over 0.5 kpc galactocentric radius bins. S c increases from 3 kpc to 4 kpc, remains 
constant between 4 kpc and 6 kpc, and decreases beyond this poi nt, from 170 pc -2 down 



to 120 M pc~ 2 at R sa \ = 7.5 kpc. It has been suggested by ( jMcKed Il989l ) that the star 
formation rate in a molecular cloud is governed by ambipolar diffusion, and subsequently by 
the surface mass density. In a molecular cloud with low surface mass density, the ionized 
fraction sustained by the ISRF penetrating the molecular cloud is high enough to prevent 
collapse and subsequent star formation due to magnetic support. In contrast, star formation 
proceeds more rapidly in molecular clouds that exhibit a high surface mass density. The de- 
cline in the average molecular cloud surface mass density with galactocentric radius beyond 
6 kpc may therefore suggest a decrease in the star formation rate with Galactocentric radius. 
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9. Discussion 



The purpose of this section is to emphasize the major assumptions and limitations of 
the results presented in this paper. Observational biases affecting the physical properties 
of molecular clouds located at different distances, as well as fundamental assumptions and 
uncertainties made in the estimation of the physical properties of molecular clouds are dis- 
cussed. 



9.1. Biases 



The top panel of Figure [9] shows the mass versus the distance to each molecular cloud. 
The mass (and radius) of the GRS molecular clouds is observed to on average increase with 
distance. This is due to i) a Malmquist bias resulting from the increasing completeness 
limit and probability of observing a massive cloud with distance, and to a minor extent ii) 
the effects of finite resolution and the method used to identify molecular clouds. Indeed, it 
is hard to define one set of parameters for CLUMPFIND (spatial and spectral smoothing, 
brightness threshold and brightness increment) that identify a molecular cloud at all emission 
levels. 



9.1.1. Modeling biases 

In order to analytically describe the effects of increasing completeness limit and probabil- 
ity of observing rarer, more massive clouds as a function of distance, we assume a power-law 
mass spectrum, 0(M) = O M -1,64 (see Section E]). The average mass of a molecular cloud 
detected by the GRS at a distance d in the GRS field is: 



grff M(j)(M)dM 



< M >(d)= j ^^:::z: m 



or 



„ t j\ M max (df- m - M min (d) - 36 

where M m j n (d) is the minimum mass that can be detected by the GRS at distance d i.e., 
the completeness limit. M min (d) is given by M min (d) = 50dl pc (see Section [6]), where dk pc is 
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the distance in kpc. M max is a physical upper limit on the masses of molecular clouds used 
to integrate Equation [181 which would diverge for M max = oo. Indeed, molecular clouds 
cannot grow infinitely due to disruptive mechanisms (e.g., turbulence, Galactic shear), the 
limited molecular scale height and limited molecular content of the Milky Way. In addition, 
the field-of-view of the GRS is limited to galactic latitudes —1° < b < 1°, such that the 
maximum mass detectable in the survey cannot exceed the limit derived in Section [6j M™ x v 
= 1.96xl0 5 <i^ pc . However, Figure [T0~| showing the maximum mass detected in each distance 
bin and the theoretical upper limit due to the limited field-of-view, demonstrates that the 
maximum mass allowed by the limited field-of-view of the GRS, M™ x v , is never reached. 
Thus, we use a constant upper bound, M max , in the integration of Equation [18j We use 
M max = 1O 6 M , the largest molecular cloud mass obser ved in the GRS, and also a typical 
high end on the mass of massive giant molecular clouds ( ISolomon et al.l 119871 ). 



9.1.2. Biases observed in the GRS 

Figure |9] shows the molecular clouds' masses versus distance. The bottom panel shows 
the trend averaged over 1 kpc distance bins. The dashed line indicates the prediction from our 
simple analytical model of the Malmquist bias. Our analytical prediction for the Malmquist 
bias is not as steep as the observed trend, and does not reproduce the observed data well. 
This could be due to several effects. First, the probability of observing rare, massive clouds 
in the solar neighborhood is low due to the small volume probed at such a low distance, 
and due to the lack of massive Galactic structure in this area. As a result the high end 
of the mass spectrum is not well randomly sampled in our data. This tends to make the 
average molecular cloud mass at low distance lower than the prediction obtained from the 
Malmquist bias. Second, the surface density of molecular gas, and thus the probability of 
observing massive molecular clouds, is not uniform in the Milky Way. 

To further investigate the cause of the observed distance/mass trend, we have performed 
a simple simulation. We have simulated a sample of molecular clouds with a mass spectrum 
4>{M) oc M -L64 that matches the Galactic surface mass density of molecular gas observed 
in Figure |HJ To that end, we have populated concentric rings of constant Galactocentric 
radius and of thickness 50 pc centered on the galactic center with molecular clouds, in order 
to match the observed Galactic surface density of molecular clouds. At each Galactocentric 
radius, the molecular cloud population in the ring is randomly sampled from a power-law 
mass spectrum of exponent —1.64, with masses ranging between 10 M & and 3xl0 6 M Q . The 
azimuth of the clouds is randomly selected from a uniform distribution between and 27r, 
and their distance and galactic longitude are computed based on their Galactocentric radius 
and azimuth. Molecular clouds that have masses greater than the completeness limit at their 
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distance (M min = 50 dl pc ) are included in the population of the ring, until the surface density 
of the ring reaches the observed surface density of molecular gas at this Galactocentric radius, 
given by Figure [HI The simulation produced a sample of 2644 molecular clouds between a 
Galactocentric radii 1.5 kpc and 10 kpc, 718 of which are located within the GRS longitude 
range. In this simulated sample of molecular clouds located inside the GRS coverage, we 
recorded the maximum and mean molecular cloud mass in each distance bin (we used 1 kpc 
distance bins), the same way we did for our GRS sample of molecular clouds. The results of 
this simulation is shown in Figure M as triangles. 

Our simple simulation of the Malmquist bias seems to reproduce the observed trends 
within the error bars. We are therefore confident that the observed trends are indeed due to a 
Malmquist bias. Effects of the molecular cloud identification algorithm may play a secondary 
role in the relation between distance and molecular cloud masses at small distances. Indeed, 
CLUMPFIND tends to detect well resolved, nearby molecular clouds as distinct objects while 
the same molecular clouds at a farther distance would be blended together and assigned to 
a single cloud. 



9.1.3. Effects of biases on molecular clouds' physical properties 

Biases mainly affect the significance of the radial profiles derived in Section [HI Because 
the masses of the molecular clouds on average increase with distance, the radial variations of 
the molecular gas content of the Milky Way seen in Figure [H] would result from biases if the 
different radial bins were dominated by molecular clouds that were on average at different 
ranges of distance. The dash-dot lines in each panel of Figure M represent < d >o.5k P c, the 
average of d in each 0.5 kpc galactocentric radius bin, versus distance. Because the Malmquist 
bias causes the average molecular cloud mass to increase with distance, < d >o.5fc pc (-Rgai) 
should give us an indication of whether or not the radial profiles shown in Figure [S] are 
dominated by biases. There is no obvious correlation between < d >o.5fc pc (-Rgai) and < M > 
(Rga) or < T, c (R ga i), indicating that the variations of < M > and S c with R ga i are genuine. 



9.2. Effects of 13 CO abundance variations 



The derivation of the GRS molecular clouds' masses is based on the assumption that 
the abundance of 13 CO relative to H2 is uniform within a molecular cloud. Owing to photo- 
dissociation and fractionation, the abundance ratio n( 13 CO)/n(H2) may however decrease 
significantly betwee n the shielded dense interiors of mole cular clouds and their diffuse, UV- 
exposed envelopes (ILisztl 120071 : iGlover &: Mac Lowl 120101 ) . Molecular hydrogen (H 2 ) being 
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self-shielded due to numerous optically thick absorption lines in the UV, it can resist the 
photo-dissociating ISRF at lower ext inctions (Ay of a few ten t h) than its CO cou nterpart, 
which can only exist for A v > 1-3 ( IGlover fc Mac Low! bold ; IWolfire et al.ll2010l ). A sig- 
nificant amount of molecular gas is therefore likely hidden in CO dark envelopes around 
molecular clouds. As a result, o u r met hod probably underestimates the masses of the GRS 
molecular clouds. IWolfire et al.l (120101 ) predicts that 30% of a molecul ar cloud's mass is in 
the form of CO-dark molecular envelope, while iGoldsmith et al.l ( 120081 ) predicts 50%. This 
is nonetheless the order of magnitude of the molecular mass not traced by CO in the Milky 
Way. 

In addition, the CO molecule tends to leave the gas phase and freeze out onto dust grains 
(thus forming a mantle) in cold, dense molecular cores. This effect could also potentially 
render some of the molecular mass inside a molecular cloud invisible to CO observations. 
Indeed, using other mole cules such as N 2 H + a nd the dust continuum as a tracer for the 
densest molecular phase, iBacmann et al.l (120021 ) have measured a CO under-abundance of 
5-15 in molecular cores of de nsity 10 5 -10 6 cm~ 3 . The under-abundance depends on the 
density (IBacmann et al.l |2002| ) . such that CO manteling on dust grain do es not occur for 



densi t ies lower than 10 4 cm 3 . The typical size of such cores is 0.05-0.2 pc (IBacmann et al. 



20001 ; iRathborne et al.ll2010l ). which at a distance of 5 kpc respresents an angular size of 
15". Thus, a core's projected area represents 1% of the GRS beam area. Therefore, we do 
not expect CO manteling on dust grain to have a significant effect on the estimation of a 
molecular cloud's mass compared to photo-dissociation effects. 



9.3. Effects of sub-thermal excitation 



In order to derive 13 CO optical depths and molecular cloud masses, the CO excitation 
temperature was assumed to be identical for 13 CO and 12 CO. Since CO is usually thermalized 
within molecular clouds due to its low dipole moment, this is a reasonable assumption. 
This assumption might however break down in the more diffuse envelopes of molecular 
clouds, where the optically thick 12 CO can remain thermalized due to radiative trapping, 
while the optically thin 13 CO is sub-thermally excited. In this case, the 13 CO excitation 
temperature would be lower than the 12 CO excitation temperature. Nonetheless, emission 
from subthermally excite d 13 CO in diffuse m olecular cloud envelopes is probably under the 
GRS detection threshold Jflever et al.ll200g[ ). 
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9.4. Uncertainties on the mass estimation 



In addition to the effects of sub-thermal excitation and abundance variations, noise 
intrinsic to the data, and errors on kinematic distances affect the accuracy of a molecular 
cloud's mass estimate. The estimation of the masses of the molecular clouds is based on the 
knowledge, at each voxel, of the excitation temperature (derived from the 12 CO brightness 
temperature), of the optical depth (derived from the excitation temperature and the 13 CO 
brightness temperature), and the distance of each molecular cloud. The 13 CO and 12 CO 
brightness temperatures are affected by gaussian noise, with a sta ndard deviation of a r mt = 
0.26 K for both the UMSB ( 12 CO) and the GRS ( 13 CO) surveys rt.Tackson et al.l l200fih. The 
accuracy of kinematic distances was discussed in Roman-Duval et al. ( 2009h . The error on 
the distance stems from the error in the velocity of a molecular cloud with respect to its LSR 
velocity, the difference being caused by cloud-to-cloud dispersion and non-circular motions 
assoc iated with spiral arms. The cloud-to-cloud dispersion amounts to 3-5 km s _1 (I Clemens 
985) and loca l velocity perturbations associated with spiral arms are of order 15 km s _1 



flClemenslll985f ). 



The error on the molecular clouds' masses due to each of those factors was estimated 
using Monte-Carlo simulations. For every molecular cloud in the GRS, we computed the 
velocity dispersion and peak brightness temperature of the 13 CO line toward each GRS 
pixel. The 13 CO line toward each pixel of the molecular cloud was then modeled by a 
gaussian line with the same peak brightness temperature and velocity dispersion. Gaussian 
noise of standard deviation cr T , = 0.26 K was added to the model molecular cloud, and its 
13 CO optical depth estimated assuming an excitation temperature of 6.3 K (the average value 
in Figure H]). We reproduced the contribution of the distance uncertainty on the mass error 
using the following method. First, a random gaussian velocity error e v of dispersion a £v = 5 
km s _1 was added to the LSR velocity of the molecular cloud, Vq, to obtain a flawed estimate 
of the molecular cloud LSR ve locity = + e v- The kinematic distance d± corresponding 
to V\ was computed using the I Clemens! (119851 ) rotation curve and the Galactic longitude of 
the molecul ar cloud, conserving the nea r /far kinematic distance assignment of the molecular 
clouds from Roman-Duval et al. ( 20091 ). The distance d\ was finally used, along with the 
13 CO optical depth of the model molecular cloud, to compute its mass. This process was 
repeated 20 times for the standard deviation of the mass of the model molecular cloud to 
converge. The error on the mass estimation of the molecular cloud was then computed as 
o M = a/ < (Mj — M c i ou d) 2 >i=o,2o, where M c i ou d is the mass of the GRS molecular cloud, and 
Mi is the mass of the i th realization of the model molecular cloud. Error bars on the mass 
estimate are specified in Table [j] and in the online complete version of the table. 
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10. Conclusion 

We have derived the physical properties (radius, mass, excitation temperature, optical 
depth, virial parameter, density, and surface density) of 580 molecular clouds identified in 
the GRS and covered by the UMSB 12 CO survey. We have derived the histograms of these 
properties, and found a power-law decrease of exponent — 1.64±0.25 for the mass spectrum 
of molecular clouds, consistent with previous results from the literature. The median virial 
parameter of nearly 0.5 suggests that molecular clouds are gravitationally bound entities. 
The range of values for the physical properties is entended to constrain numerial models of 
molecular cloud formation and evolution. 

We have found a tight power-law correlation of exponent 2.36±0.04 between the radii 
and masses of molecular clouds. Based on this correlation, we have deduced that the fractal 
dimension of the ISM must be of order 2.36 in the molecular phase, and in the range of 
spatial scales covered by molecular clouds. The correlation between molecular cloud radii 
and masses also allowed us to derive masses for an additional 170 GRS molecular clouds not 
covered by the UMSB 12 CO survey, for which excitation temperatures and optical depths 
could not be derived. Based on the 750 molecular clouds, we have examined the Galactic 
surface mass density of molecular clouds in the Galactic region covered by the GRS. The 
azimuthally averaged Galactic surface mass density of molecular gas is enhanced by a factor 
of 5 between galactocentric radii of 4 and 5 kpc. This supports previous observations that 
most of the molecular cont ent of the Milky Way is contained in a ring located 4-5 kpc away 



from the Galactic center ( jClemens et al.lll988[ ). In addition, the two-dimensional map of 



the Galactic surface mass density of molecular gas is enhanced along the positions of the 
Scutum-Crux and Sagittarius arms (also detected with other tracers such as star counts, H II 
regions, magnetic fields, electron density, etc ...). We have found no enhancement of molecu- 
lar gas along the Perseus arm, which could be due to uncertainties in the kinematic distance, 
near/far blending, or the lack of completeness at this distance (> 10 kpc). Nonetheless, the 
molecular gas enhancements observed at the assumed positions of the Scutum and Sagittar- 
ius arms may suggest that molecular clouds form in spiral arms, and are disrupted in the 
inter-arm space. 

Last, we have found that the CO excitation temperature of molecular clouds decreases 
with galactocentric radius, suggesting a decrease in the star formation rate away from the 
Galactic center. The excitation temperature of molecular clouds is also marginally enhanced 
at a galactocentric radius of 6 kpc, which, in the GRS longitude range, corresponds to the 
inferred position of the Sagittarius arm. This marginal increase may be related to star for- 
mation activity in the Sagittarius arm. 
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Fig. 1. — Correlation between the masses and radii of molecular clouds. 
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Fig. 2. — Positions of the GRS molecular clouds covered by the UMSB 12 CO survey (black 
crosses), and outside of the UMSB coverage (red crosses). The dashed circles indicate galac- 
tocentric radii R, ga i = 1-8 kpc, by steps of 1 kpc. The solid circle indicates the solar circle. 
The white area corresponds to the portion of the Galactic plane covered by the GRS. The 
artefact produced by the alignment of clouds in an arc of circle is due to the tangent point. 
Indeed, molecular clouds with radial velocities greater than the radial velocity of the tangent 
point (due to uncertainties in the rotation curve and non-circular motions) were assigned 
the distance of the tangent point. 
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Fig. 3. — Histograms of the physical properties of molecular clouds. In the top two panels, 
the solid line indicates the best fit to the radius and mass spectra. 
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Fig. 4. — Histograms of the excitation temperature (top) and 13 CO optical depth (bottom) 
of the molecular clouds averaged over voxels where the brightness temperature is greater 
than 4(7 = 1 K. 
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Fig. 5. — Galactic surface mass density of molecular clouds. The contour levels are 10, 20, 
30, 40, 50, 60, 70, 80, 90% of the maximum density (5.3 x 10 6 ). The dashed lines indicate 
galactocentric radii of 1, 2, 3, 4, 5, 6, 7, and 8 kpc. The solid line indicates the solar 
circle. The white area corresponds to the portion of the Galactic plane covered by the GRS. 
The green, yellow, red, and blue lines represent the 3 kpc arm, t he Scu t um- C rux arm, the 
Sagittarius arm, and the Perseus arm from the four-arm model by I Vailed (119951 ) respectively. 



Table 1. Catalog of molecular clouds's masses and physical properties. 

GRS molecular cloud 1 b V LS R AV R M 6 M n(H 2 ) T cx t( 13 CO) 

(° ) (° ) (kins" 1 ) (kms- 1 ) (pc) (M ) (M ) (c m - 3 ) (K) 



GRSMC 


G053 


59+00 


04 


53 


59 


0.04 


23 


74 


1 


86 


17.5 


2 


35 


X 


10 s 


0.566 X 


10 5 


161.0 


8 


28 


1 


78 


244.4 


0.04 


GRSMC 


G049 


49-00 


41 


49 


49 


-0.41 


56 


90 


9 


12 


9.8 


1 


81 


X 


10 5 


0.445 X 


10 5 


707.7 


9 


97 


1 


47 


601.6 


0.73 


GRSMC 


G018 


89-00 


51 


18 


89 


-0.51 


65 


82 


2 


80 


12.4 


1 


41 


X 


10 5 


0.457 X 


10 5 


271.9 


9 


72 





89 


292.4 


0.11 


GRSMC 


G030 


49-00 


36 


30 


49 


-0.36 


12 


26 


4 


56 


1.7 


7 


82 


X 


10 s 


3.72 X 


1() 2 


617.7 


5 


97 


1 


82 


89.3 


7.25 


GRSMC 


G035 


14-00 


76 


35 


14 


-0.76 


35 


22 


4 


89 


5.4 


4 


95 


X 


10 s 


1.21 X 


10 4 


1175.8 


8 


69 


2 


21 


548.0 


0.42 


GRSMC 


G034 


24+00 


14 


34 


24 


0.14 


57 


75 


5 


66 


12.3 


1 


44 


X 


10 s 


0.406 X 


10 5 


281.0 


6 


48 


1 


83 


300.6 


0.45 


GRSMC 


G019 


94-00 


81 


10 


94 


-0.81 


42 


87 


2 


62 


8.2 


5 


54 


X 


10 s 


1.64 X 


10 4 


366.9 


8 


90 





86 


261.8 


0.17 


GRSMC 


G038 


94-00 


46 


38 


94 


-0.46 


41 


59 


2 


78 


22.8 


4 


88 


X 


10 s 


1.03 X 


10 5 


150.2 


8 


27 


1 


77 


297.8 


0.06 


GRSMC 


G053 


14+00 


04 


53 


14 


0.04 


22 


04 


2 


26 


4.0 


1 


09 


X 


10 s 


0.397 X 


10 4 


651.0 


7 


05 


2 


00 


222.9 


0.30 


GRSMC 


G022 


44+00 


34 


22 


44 


0.34 


84 


52 


2 


62 


7.5 




22 


X 


10 5 


1.20 X 


10 4 


448.2 


7 


99 


1 


29 


293.3 


0.16 


GRSMC 


G049 


39-00 


26 


49 


39 


-0.26 


50 


94 


3 


51 


12.3 


1 


57 


X 


10 s 


0.409 X 


10 5 


309.6 


9 


17 


1 


56 


330.2 


0.16 


GRSMC 


G019 


39-00 


01 


19 


39 


-0.01 


26 


72 


3 


63 


6.3 


4 


93 


X 


10 s 


1.85 X 


10 4 


732.0 


7 


11 


2 


05 


399.1 


0.27 


GRSMC 


G034 


74-00 


fifi 


34 


74 


-0.66 


46 


69 


4 


09 


8.3 


9 


77 


X 


10 5 


2.72 X 


10 4 


629.5 


8 


30 


1 


2(1 


452.3 


0.23 


GRSMC 


G023 


04-00 


41 


23 


04 


-0.41 


74 


32 


4 


14 


15.8 


2 


23 


X 


10 s 


0.679 X 


10 5 


208.9 


9 


16 





92 


285.6 


0.20 


GRSMC 


G018 


69-00 


06 


18 


(19 


-0.06 


45 


42 


3 


85 


21.8 


4 


78 


X 


10 s 


1.24 X 


10 5 


169.6 


7 


78 


1 


52 


320.7 


0.11 


GRSMC 


G023 


24-00 


36 


23 


24 


-0.36 


77 


30 


2 


65 


15.5 


1 


95 


X 


10 5 


0.683 X 


10 5 


190.9 


8 


64 


1 


01 


257.3 


0.09 


GRSMC 


G019 


89-00 


56 


19 


89 


-0.56 


44 


14 


3 


68 


8.9 


8 


63 


X 


10 5 


2.82 X 


10 4 


446.3 


7 


83 


1 


21 


345.3 


0.23 


GRSMC 


G022 


04+00 


1!) 


22 


04 


0.19 


50 


94 


5 


53 


10.8 


8 


48 


X 


10 s 


1.93 X 


10 4 


248.3 


7 


17 


1 


15 


232.3 


0.63 


GRSMC 


G018 


89-00 


66 


18 


89 


-0.66 


64 


12 


3 


66 


12.6 


1 


59 


X 


10 5 


0.404 X 


10 5 


294.5 


8 


90 





97 


320.8 


0.17 


GRSMC 


G023 


34-00 


21 


23 


34 


-0.21 


81 


12 


4 


09 


13.1 


1 


63 


X 


10 5 


0.509 X 


10 5 


267.4 


8 


87 





97 


303.6 


0.22 


GRSMC 


G034 


99+00 


34 


34 


99 


0.34 


53 


07 


3 


57 


11.8 


1 


28 


X 


10 s 


0.352 X 


10 5 


289.6 


7 


35 


1 


41 


295.0 


0.19 


GRSMC 


G029 


64-00 


(il 


29 


64 


-0.61 


75 


60 


3 


79 


13.1 


9 


58 


X 


10 s 


3.01 X 


10 4 


156.9 


6 


61 


1 


37 


178.2 


0.32 


GRSMC 


G018 


94-00 


26 


18 


94 


-0.26 


64 


55 


2 


72 


10.4 


1 


17 


X 


10 5 


0.402 X 


10 5 


381.0 


8 


79 


1 


27 


344.2 


0.11 


GRSMC 


G024 


94-00 


l(i 


24 


94 


-0.16 


47 


12 


4 


40 


9.5 


9 


21 


X 


10 5 


2.14 X 


10 4 


398.5 


8 


57 


1 


00 


327.5 


0.32 


GRSMC 


G025 


19-00 


26 


25 


19 


-0.26 


63 


70 


2 


30 


5.0 


1 


73 


X 


10 s 


0.438 X 


10 4 


501.7 


7 


90 


1 


48 


218.4 


0.25 



-38 - 




Fig. 6. — Excitation temperature of molecular clouds versus galactocentric radius. The red 
crosses represent the maximum excitation temperature of each molecular cloud. The black 
crosses represent the mean excitation temperature of each molecular cloud, averaged over 
voxels with a brightness temperature greater than 4a = 1 K. In the lower panel, the excitation 
temperature of GRS molecular clouds were averaged over 0.3 kpc galactocentric radius bins. 
The CO excitation temperature can be lower than the gas temperature if the local number 
density is lower than the critical density. The maximum excitation temperature in the clouds 
shoud be close to the gas temperature. The dashed blue and green lines represent the gas 
temperature predicted from thermal equilibrium between cooling by CO rotational emission 
and heating by Galactic cosmic rays (blue) and photo-electric effects from dust grains (green) 
respectively. Enhancements of 50% in the GCR flux and the ISRF were included in the blue 
and green curves respectively at K ga i = 6.4 kpc to explain the observed enhancement in the 
CO excitation temperature. 
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Fig. 7. — Distributions of excitation temperatures of "active" (i.e., containing an Hll region) 
and "quiescent" molecular clouds in red and blue respectively. On average, active molecular 
clouds have slighlty higher temperatures. 
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Fig. 8. — Radial variations of i) the Galactic surface mass density of molecular clouds (left), 
ii) of the molecular cloud masses (middle), and iii) of the molecular clouds' surface mass 
density (right). The dashed line represents the average < d >o.5k P c in each 0.5 kpc bin. 
< d > . 5fc pc is an indicator of how much the average physical properties of molecular clouds 
within each radial bin are affected by biases (see section I9~TT) . 
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Fig. 9. — Molecular clouds' individual masses and masses averaged within 1 kpc distance 
bins versus distance. The dashed lines represent the expected trend from the analytical 
description of the Malmquist bias predicted for the GRS. The triangles are predictions from 
a numerical simulation of the Malmquist bias in the GRS. 
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Fig. 10. — Maximum molecular cloud mass detected in each 1 kpc distance bin as a function 
of distance. The dashed line represents the maximum mass allowed by the field-of-view of 
the GRS. 



